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We determine elastic and coupled-channel amplitudes for isospin-1 meson-meson scattering in 
P-wave, by calculating correlation functions using lattice QCD with light quark masses such that 
= 236 MeV in a cubic volume of ~ (4fm)®. Variational analyses of large matrices of correlation 
functions computed using operator constructions resembling tttt, KK and qq, in several moving 
frames and several lattice irreducible representations, leads to discrete energy spectra from which 
scattering amplitudes are extracted. In the elastic tttt scattering region we obtain a detailed energy- 
dependence for the phase-shift, corresponding to a p resonance, and we extend the analysis into the 
coupled-channel KK region for the first time, finding a small coupling between the channels. 


The study of hadron spectroscopy from first principles 
QCD is entering a new stage of development where the 
relationship between the discrete spectrum of the theory 
in a finite-volume and the infinite-volume scattering am¬ 
plitudes is being practically utilized to study resonances. 
The tool which allows us access to the spectrum is lattice 
QCD in which the quark and gluon fields are considered 
on a finite-grid of points with only systematically improv¬ 
able approximations being made. 

Predictably it is the simplest resonant scattering chan¬ 
nel which has attracted the greatest initial interest [1-5] , 
that of TTTT with isospin=l, in which a low-lying elas¬ 
tic vector resonance called the p appears. These works 
have made use of the formalism relating the discrete spec¬ 
trum at rest and in moving frames to elastic scattering 
amplitudes, which has been in place for many years [6- 
10]. Recently the extension to coupled-channels has been 
presented [11-14], and the first lattice QCD study of a 
coupled-channel system, that of nK, pK in S, P and D- 
waves, has appeared [15, 16], showing that the energy 
dependence and resonant content of the scattering ma¬ 
trix for such a system can be extracted from finite volume 
spectra. 

To date virtually all determinations of hadron scatter¬ 
ing amplitudes in lattice QCD calculations have worked 
with artificially heavy u, d quark mass values, a choice 
which leads to heavier than physical pseudoscalar mesons 
- this reduces the computational cost, allowing calcula¬ 
tions in smaller volumes (where remains large), and 
pushes up in energy the thresholds for multihadron scat¬ 
tering such as TTTTTTTr, for which a finite-volume formalism 
is not yet in place (but see Refs. [17-20] for progress). 
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The Hadron Spectrum Collaboration previously com¬ 
puted TTTT scattering using 391 MeV pions [21, 22], ex¬ 
tracting detailed spectra of QCD eigenstates from vari¬ 
ational analysis of two-point correlation functions com¬ 
puted in several moving frames in three different vol¬ 
umes. By obtaining a significant number of energy levels 
in the elastic scattering region they were able to map out 
the energy dependence of the scattering amplitude and 
show that there is a narrow p resonance barely above tttt 
threshold. 

In this paper we deliver an extension of the work pre¬ 
sented in [22], utilizing a smaller u,d quark mass, corre¬ 
sponding to a pion mass of 236 MeV, in a large box of 
spatial extent ^ 3.8 fm. Going beyond what was done 
before, we also consider the effect of including KK-like 
operators into the variational operator basis - the en¬ 
larged basis allows us to determine energy levels above 
the KK threshold, and to extract first estimates within 
QCD of the coupled-channel tttt, KK scattering matrix 
with 1=1, = 1~. 


I. CALCULATING THE FINITE-VOLUME 
SPECTRUM 

The results to be presented in this paper come 
from a calculation using a single ensemble of 
anisotropic Clover gauge-field configurations of volume 
(L/os)^ X (T/at) = 32^ x 256, with spatial lattice spac¬ 
ing Us 0.12fm, and temporal lattice spacing a* * * § = Us/^ 
with ^ ~ 3.5. The 2 -|- 1 flavors of dynamical quarks 
have strange quark mass tuned to approximate the phys¬ 
ical strange quark [23, 24] and degenerate u, d quarks 
with mass parameter atirii = —0.0860 corresponding to 
a pion mass ~ 236 MeV [25-27]. The large volume and 
time extent, rriTrL ^ 4.3, and ^ 10, ensure that ex¬ 
ponentially suppressed polarization and thermal effects 
will be negligible. Correlation functions are computed 
on 469 configurations, typically utilizing multiple time 
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sources on each to increase statistics. 

Our approach is to determine the spectrum from a ma¬ 
trix of two-point correlations functions constructed using 
a basis of hadronic operators at source and sink. Our 
basis, which is described extensively in previous publica¬ 
tions [21, 28-30], contains both “single-meson-like” oper¬ 
ators of the form '0r D ... Dip, and “meson-meson-like” 
operators of the form O'l'(pi) Ol'(p 2 )- 

(p) is a variationally optimized combination of “single¬ 
meson-like” operators, capable of interpolating a sta¬ 
ble pseudoscalar with momentum p. The use of vari¬ 
ous meson momenta, pi,P 2 , at a fixed total momentum, 
P = Pi -\-p 2 , furnishes the “meson-meson-like” part of the 
operator basis, where the generalized Clebsch-Gordan co¬ 
efficient, C{pi,p 2 ) ensures the operator has the desired 
properties under rotations and parity. 

In order to efficiently compute the correlation func¬ 
tions, whose Wick contractions include quark-antiquark 
annihilation on a timeslice, we make use of distilla¬ 
tion [31]. The distillation smearing operator is con¬ 
structed as an outer product of the lowest 384 eigen¬ 
vectors of the gauge-covariant laplacian on each times¬ 
lice. The light and strange quark propagators required 
to evaluate the correlation functions are the solutions of 
the Clover-Dirac equation using these eigenvectors as a 
source^. To our knowledge, this is the largest lattice 
volume on which the full distillation method has been 
applied. 

We form the possible combinations of creation and an¬ 
nihilation operators at source and sink to construct a 
matrix of correlation functions. This can be analyzed 
variationally [36-38] by solving a generalized eigenvalue 
problem, C{t)v'' = X„{t)C{to)v'', where the eigenval¬ 
ues An(t) give information about the spectrum, and the 
eigenvectors provide the optimal linear combination of 
basis operators to interpolate state |n). Details of our 
implementation may be found in Refs. [28, 29]. 

Through calculating the discrete spectrum in moving 
frames we may better constrain our description of scat¬ 
tering amplitudes - each moving frame gives the bound¬ 
ary a different geometry, which lead to modified quan¬ 
tization conditions relating amplitudes to the spectrum. 
In a finite, spatially periodic, cubic volume, there is a re¬ 
duced symmetry with respect to that of an infinite space. 
For systems with no overall momentum, parity is a good 
quantum number and a double cover of the octahedral 
group, Oj^, describes the symmetries of the system. The 
infinite volume partial waves are subduced into various 
irreducible representations, or “irreps”, of the octahedral 
group. Systems may be considered having non-zero over¬ 
all momentum, P = ^ [n^,, n^, n^], which satisfy peri¬ 
odic boundary conditions if nx,ny,n 2 are integers. Such 


^ The large number of propagators are very efficiently computed 
using an Adaptive Multi-Grid solver [32, 33] for the light quarks 
on CPUs, and the strange quark propagators are efficiently com¬ 
puted using Graphical Processing Units [34, 35]. 



FIG. 1. Pseudoscalar tt, K and rj meson dispersion relations 
obtained on this lattice for jPp (^) < 4 as described in the 
text. The statistical uncertainties on the energies are smaller 
than the symbols. The curves correspond to dispersion re¬ 
lation fits of the form (atE)^ = -|- ^ nP with 

parameters presented in Table I. 

systems have a further reduced symmetry, relative to the 
rest-frame, and are described by the little groups, LG(.P), 
which are the subgroups of the octahedral group whose 
transformations leave P unchanged. The irreps of LG(P) 
typically have an enlarged angular momentum content 
with respect to the system at rest, and parity is not a 
good quantum number. These concepts are discussed in 
detail in Ref. [30] with applications relevant to tttt scat¬ 
tering further developed in Refs. [21, 22]. 

The masses (and energies at nonzero momentum) 
of the stable pseudoscalars, tt, K and rj, are ob¬ 
tained from the lightest eigenstates in the irreps 
[000] and |P| > 0, A 2 , and are presented in Fig¬ 
ure 1, along with dispersion relation fits of the form 

(atE)'^ = {atmY -h {^ ^2]a ) where ^ is allowed to 
vary for each pseudoscalar species - Table I shows the 
masses and extracted ^ values. We observe a reasonable 
level of consistency in ^ between species, and we choose 
to use the value from the pion, ^ = 3.453(6), in the rest of 
this paper. The variational analysis that yields these en¬ 
ergies also provides the eigenvectors used when forming 
the optimized operators D that are used in the “meson- 
meson-like” constructions. 

Utilizing these methods, we computed 1 = 1, G-parity 
positive spectra for all irreps^ containing = 1“ with 
|p|2 < 4 (^) . The operator basis used to construct 
the correlation matrices for each lattice irrep. A, is pre¬ 
sented in Table II, where we see that the dimension of 
the correlator matrices to be considered is never smaller 
than 10 and may be as large as 37. As an example of 


^ A table of the subductions of tttt (and equivalently KK) partial- 
waves into these irreps appears as Table III in [22]. 
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threshold 

atE 


at m 


xVlVdof 

TTTT 

0.07856(26) 

TT 

0.03928(18) 

3.453(6) 

1.4 

TrTTTrTr 

0.15712(52) 

K 

0.08344(7) 

3.462(4) 

1.4 

KK 

0.16688(14) 

V 

0.09299(56) 

3.468(20) 

0.73 

rjTvn 

0.17155(62) 





nKK 

0.20616(19) 


TABLE I. The masses of stable particles and anisotropy, 
obtained from dispersion fits using five momenta, as described 
in the text, and presented in Figure 1. Threshold energies 


relevant to 
presented. 

I = 1, J’" 

= 1 meson 

-meson scattering also 

[000] Tf 

[100] Ai 

[110] Ai 

[111] Ai 

[200] Ai 

TTOOl TTOO-l 

TTOll TTO-l-l 

TTlll 7r_i_i_i 

TTOQQ TTlOO 

TTO-IO TTllO 

TTO-l-l TTlli 

TT-lOO 7r200 

Tl'OOO Tl'llO 

TTOO-I TTlll 

TT-liO 7r200 

TTOOO TTlll 

TTlOO TTOll 

TT-lll 7r200 

TTOOO 7r200 

A'ooi Aoo-i 

A’ooo Aioo 
Ao-io Alio 

Aooo Alio 
Aoo-i Am 

o o 

o o 


tpTtp X 26 

ipFip X 10 

ipFip X 13 

iljFij; X 21 

tpTip X 18 



[100] Aa 

[110] Bi 

[110] B2 

[111] Aa 

[200] Aa 

TTO-IO TTllO 

TI'O-l-l TTlll 

Tl'OlO TTIOO 

TI'Ol-l TTlOl 

TT-ilO 7r200 

TTOO-l TTiii 

TI'Ol-l TTlOl 

TTlOO TToil 

Tl'-lll 7r200 

TTl-lO TTiio 

TI'l-l-l TTlll 

ijjF'ip X 29 

tpFip X 29 

ipFtp X 29 

tpTip X 35 

ipFip X 29 



[100] Bi 


[100] Aa 



TTO-lO TTllO 


TTO-l-l TTill 



tjjFij} X 9 


ipF'ip X 9 



TABLE 11. The operator bases used in each lattice irrep in this 
calculation. For each irrep we list the “TTTr-like” and “KK- 
like” operators (Appendix B contains details of the KK op¬ 
erator construction) that were used as well as the number of 
“single-meson-like” operators. We use a notation which indi¬ 
cates the momentum (in units of 27r/L) of the pseudoscalar 
meson operators, recalling that the directions of momentum 
are summed over with generalized Clebsch-Gordan weights to 
ensure the operator lies in the stated irrep [21, 22]. 


the result of variational analysis, in Fig. 2 we present the 
first five An(t) from the A^ = irrep which illustrates 
the quality of the spectra obtained. Higher energy levels 
are extracted, but are not used in the scattering analysis 
which follows. 

The spectra extracted are shown in Figures 3, 4 and 
5, where we display the energy levels transformed to 
their cm-frame values (points with statistical error bars), 
the scattering thresholds (short-dashed horizontal lines) 
and the spectrum of non-interacting tttt and KK levels 
(curves). As can be seen in Table II, we have not included 
any operators featuring pseudoscalar operators with mo¬ 
mentum, IpI^ >4. The long dashed curves indicate those 


non-interacting levels for which we have not included the 
corresponding operator construction. The high-lying ex¬ 
tracted levels displayed by gray circles lie in an energy 
region in which we have not included sufficient operators 
to reliably extract the entire spectrum, and these levels 
will not be used in the analysis which follows. While the 
Ai little-group irreps contain no subductions from the 
= I'*' partial-wave, the E and B irreps do, and this is 
likely to be the origin of the consistent “additional” level, 
present near atEcm ~ 0.22 in each pane of Figure 4, be¬ 
ing due to a positive parity bi resonance. These states 
are observed to have large overlap onto “single-meson” 
operators subduced from = 1+. 

The light gray dashed horizontal lines in Figures 3, 
4 and 5 show the multihadron thresholds, 47r, ijirn and 
ttKK. Note that we have not included operators resem¬ 
bling these in our variational basis, nor have we plotted 
the corresponding non-interacting levels in the Figures - 
these will lie at a higher energy than the threshold. Ex¬ 
perimentally these channels have very small amplitudes 
in e“'"e“ annihilation until several hundred MeV above 
threshold [39, 40], so we do not expect them to play a 
significant role - we will discuss this in greater detail 
later in this paper. 

Before attempting to determine meson-meson scatter¬ 
ing amplitudes from the spectra presented in Figures 3, 
4 and 5, we will present a brief illustration of the im¬ 
portance of using a sufficiently diverse basis of opera¬ 
tors in variational analysis. Figure 6 shows the spec¬ 
trum extracted in the [000] Tj“ irrep using five different 
choices of operator basis. The histograms show the rel¬ 
ative strength of overlap (nlC^lO) for the various oper¬ 
ators in the basis^. The leftmost column is our largest 
basis, the one presented in Table II which contains three 
TTTr-like operators, one KK-\ike operator, and 26 “single¬ 
meson-like” operators, 19 of which are subduced from 
= 1“ constructions, 6 from = 3“ subductions 
and one from a = 4“ subduction. The second col¬ 
umn lacks the KK-like operator and is seen to give a 
reasonably consistent spectrum with the exception of the 
level which had large overlap onto the KK-\ike operator. 
The third column uses only the “single-meson-like” op¬ 
erators, lacking any TTTr-like or KK-\ike constructions - 
the only low lying state extracted appears to be some 
sort of crude average of the two lowest lying states. The 
fourth and fifth columns, which exclude “single-hadron¬ 
like” operators also provide poor determinations of the 
spectrum. It appears, as one might expect for a system 
in which we expect a narrow resonance, usually thought 
of as a tightly bound qq state, strongly coupled to tttt, 
an accurate spectrum cannot be determined without in¬ 
cluding both “single-meson-like” operators and TTTT-like 
operators. A simple argument explaining these observa¬ 
tions (illustrated using an in-flight irrep) was previously 


® see [28, 29] for further details of the normalization of such over¬ 
laps 
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-V atE = 0.12488(40) 
xVNioS = 1.38 


10 15 20 25 t/at 

atE = 0.20275(37) 


Ao 



X /-^dof = 0.87 1.20 


atE - 0.21384(65) 
^4 X^/N^of - 1.20 
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FIG. 2. Principal correlators, \n{t), from the variational analysis of a 30 x 30 matrix of P = [000], correlation 

functions, with to = 6at. The lowest five states are plotted as ^n{t), along with a two-exponential fit determining the 

energy, E„. The bottom right panel shows the corresponding effective masses, atmeff = ^ log ^ "'ith St = Sat- 



FIG. 3. Finite volume spectra obtained in the T/~ and moving frame Ai irreps. Short dashed lines denote meson-meson and 
multimeson thresholds, with tttt in red and KK in green. Solid curves are allowed non-interacting meson-meson energies in 
finite volume, corresponding to operator constructions that were included in our basis, while the long-dashed curves are those 
that were not included in our basis. The points show the energy levels with their statistical errors as extracted from the lattice 
QGD correlation functions, with those in black being the ones used in the amplitude analysis to follow and those in grey not 
used, as described in the text. 


given in [22]. In Appendix A we discuss the result of 
performing a phase-shift extraction using the spectrum 
extracted using only “single-meson-like” operators. 

We briefly comment that in the [000] irrep, there 
are extracted levels near a* Acm ~ 0.28 which have signifi¬ 
cant overlap onto the {jirtp operators, these likely indicate 
the mass scale of the higher excited vector resonances. 
We also find levels that have significant overlaps with op¬ 
erators subduced from continuum = 3“ and = 4“ 


above atEcm = 0.33, suggesting p^, p 4 resonances. 


II. DETERMINING SCATTERING 
AMPLITUDES 

The relationship between a two-body coupled-channel 
scattering t-matrix and the discrete spectrum for an irrep 
P, A in a finite, periodic, L x L x L volume is provided 
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FIG. 4. As Fig. 3 for the E- and B-type little-group irreps which have — 1 


as the lowest subduced partial wave. 




FIG. 5. As Fig. 3 for the B-type little group irreps with 
P = [100] which have = 3“ as the lowest subduced partial 
wave. 


by the equation, 


det 
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nn' 


tw = i 

sindf, and scattering can be described by a 
single real function, the scattering phase-shift, 5^{Ecm)- 
For a given t-matrix, the solutions of Eq. 1 provide 
the finite volume spectrum, {i?n}, in each lattice ir- 
rep A with some overall momentum P. In the elastic 
case, if higher partial-waves have negligibly small am¬ 
plitudes, as one expects at low energies, the equation 
can be solved for 6i{En) for each calculated E„. In a 
two-channel scattering problem there are three unknown 
functions of energy to determine for each partial-wave^ 
so for a given level E^ this equation is underconstrained. 
If higher partial waves are not negligible, there will be 
still further unknowns. Fortunately, we are able to ex¬ 
tract multiple energy levels in many irreps and these can 
be simultaneously used to constrain the scattering am¬ 
plitude as a function of energy. By parameterizing the 
energy-dependence of the t-matrix, we can minimise a 
function describing the difference between the calcu¬ 
lated spectrum and the spectrum given by solutions of 
Eq. 1 for the t-matrix parameterization®. This method 
was first applied to a coupled-channel situation using lat¬ 
tice QCD energy levels in Ref. [15] and further details of 
this method and our implementation may be found in 
Ref. [16]. 


' iPi{Ecm) {Ecm) {Su'Snn' + )) 


= 0 . 
'( 1 ) 


A. Elastic TTTT scattering 


where the determinant is over the channel indices 
i and the partial-waves, £, subduced into irrep A. 
Pi{Ecm) = 2fci/i?cm is the phase space for channel i, 

and the finite-volume dependent matrix with 

qi = kiL/2'K where ki is the cm momentum in channel f, 
is described in Eq. 7 of Ref. [16] and Eq. 89 of Ref. [7]. 
This expression was derived in Refs. [11-14], and in the 
case of a single open channel, reduces to the conditions 
presented earlier in [6] and [7, 8]. In the elastic case. 


We first study the elastic region, by considering only 
those levels extracted below the 47r threshold, which lies 
slightly below the KK threshold. We will initially pro¬ 
ceed assuming that only tttt scattering in P-wave is rele- 


^ three independent parameters are required to describe a unitary, 
time-reversal invariant, two-channel t-matrix 
® The explicit form of the provided in Eq. 9 of ref. [22] 
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FIG. 6. The spectra obtained from various choices of operator basis in the lattice irrep. The leftmost column contains 
all of the operators we considered, including “single-meson-like” operators subduced from J = l(orange), J = 3(blue) and 
J = 4(grey). The remaining columns use fewer operators as indicated. The histograms show the suitably normalized magnitudes 
of the contributions of each operator to each energy level, {n|C)^|0), obtained from the variational solutions. See [28, 29] for 
more details. 


vant in this energy region, and later show that indeed the 
TTTT F-wave amplitude and higher play no significant role. 
When partial waves above £ = 1 are negligible, then us¬ 
ing Eq. 1 one can obtain an estimate of <5i(i?cm) at each 
calculated value of as plotted in Figures 3 and 4. 
These phase-shift values are plotted in Figure 7, where 
we see a clear resonant behavior above tttt threshold. 

In order to describe the resonant content of the scatter¬ 
ing amplitude we may explore energy-dependent param- 
eterizations. We will consider various choices of energy- 
dependent parameterization in the minimization de¬ 
scribed above and will later discuss their pole content, 
finding that all choices capable of describing the finite 
volume spectrum have a pole at the same location in 
the complex-energy plane, corresponding to a single res¬ 
onance. 

In elastic scattering, the Breit-Wigner parameteriza¬ 
tion is commonly used to describe isolated resonances - 
in our case, with s = i?cmi takes the form 


t(s) 


1 v^r(s) 

p{s) — s — i^/sT{s) ’ 


( 2 ) 


with the energy dependent width, r(s) = includ¬ 

ing a coupling constant, gu, and the threshold behavior 
required in P-wave scattering. Attempting to describe 
22 levels below 47r threshold, we find the following pa¬ 
rameters. 


1 

xV^dof = ti = 1.25 , 
where the first set of errors describes the statistical un¬ 
certainty and the second comes from varying the pion 
mass and anisotropy, within their uncertainties. The 
matrix illustrates the statistical correlation between pa¬ 
rameters, which in this case is seen to be very small. The 
corresponding (5i(Fcm) is plotted in Figure 8. 

Modifications to the Breit-Wigner form which tame 
the barrier behavior at higher energies can be con¬ 
sidered [22, 41] - fits to 22 levels with these forms yield 
barely improved values and values of mu and gn that 
are statistically compatible with those given above. Re¬ 
stricting the energy region being described by the Breit- 
Wigner of Eq. 2 to 0.117 < atE^m < 0.146, i.e. excluding 
the tails of the resonance, leaves 14 energy levels - fitting 
these also leads to mji,gji values compatible with those 
given above. The corresponding phase-shifts for these 
modified fits are plotted in Figure 9. 

A more flexible parameterization scheme is provided 
by the iX-matrix, which automatically satisfies unitar- 
ity in the single-channel and coupled-channel cases. Our 
implementation is presented in [16] and reads, for Awave 
scattering, 

where Kij(s) is a real function, and Iij{s) is the 
Chew-Mandelstam function whose imaginary part above 


mR = 0.13171 (36) (6) • Ot" 
gR = 5.691 (70) (25) 


0.04 

1 

( 3 ) 
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FIG. 7. P-wave tttt isospin-1 elastic scattering phase-shift extracted from energy levels plotted in Figures 3 and 4, assuming 
P-wave and higher partial-wave amplitudes are negligible in this energy region. Two points whose phase-shift values have 
rather large error bars are not shown. Grey dashed vertical line shows the ttti'R'it threshold. 



FIG. 8. P-wave tttt isospin-1 elastic scattering phase-shift. 
Points as in Figure 7. Gurve shows the Breit-Wigner descrip¬ 
tion whose parameters are given in Eq. 3. 

thresholds, lm/y (s) = —5ij Pi{s), ensures unitarity, and 
whose real part is defined by a dispersive integral that 
ensures that t{s) has no pseudothreshold branch point 
(at s = 0). 

In single-channel tttt scattering with £ = 1, the K- 
matrix is just a single function, and a convenient form 
is 

2 ^ / \ ^ 

Kis) = —§■ - 

with So = 4mJ. Along with a suitable subtraction in the 
dispersive integral for J(s) so that Re J(s = m^) = 0, this 
gives a t{s) behavior that is similar to a Breit-Wigner, but 
with the polynomial allowing more freedom in the en¬ 
ergy dependence. The 22 energy levels below 47r thresh¬ 



FIG. 9. P-wave tttt isospin-1 elastic scattering phase-shift. 
Points as in Figure 7. Gurves show the Breit-Wigner fit to 
the whole elastic region (grey), a Breit-Wigner with Hippel- 
Quigg [41] barrier corrections (orange), and a Breit-Wigner 
description of a narrower energy region around the resonance 
peak (red). 

old have been described by this form for three choices, 
N = 0,1,2, and without any polynomial term at all - 
the results are presented in Table III. There is negli¬ 
gible improvement in the x^/-^dof adding terms linear 
or quadratic in s, and the corresponding parameters are 
found to possess an increasingly large degree of correla¬ 
tion. The phase-shifts corresponding to the fits in Ta¬ 
ble III are plotted in Figure 10. 

The assumption that tttt P-wave scattering plays a 
negligible role in determining the spectrum in the elas¬ 
tic region, as was assumed in the previous analysis, can 
be tested using the energy levels we have determined. 
The irreps [100] Bi and i ?2 have = 3“ as their 














N 

atm 

9 

70 /at 

71 /a? 

72 /at 

xVA'dof 

- 

0.13172(36) 

0.4475(52) 

- 

- 

- 

27.0/(22 - 2) = 1.3 

0 

0.13164(36) 

0.4611(66) 

5.4(17) 

- 

- 

16.8/(22-3) = 0.88 

1 

0.13161(37) 

0.4677(82) 

-3.3(67) 

2.6(22) 

- 

15.6/(22-4) = 0.86 

2 

0.13165(37) 

0.4679(89) 

-21.5(74) 

16.6(24) 

-2.4(4) 

14.8/(22-5) = 0.87 


TABLE III. if-matrix descriptions of the elastic spectrum using Eq. 5. 



FIG. 10. P-wave tttt isospin-1 elastic scattering phase-shift. 
Points as in Figure 7. Curves illustrate the if-matrix descrip¬ 
tions given in Table III, plus a TF-matrix description using 
the simple phase-space, p, in place of the Chew-Mandelstam 
function. 


leading partial-wave, and we have determined one low- 
lying energy level in each, as shown in Fig. 5. Us¬ 
ing Eq. 1 to obtain the phase-shift, under the reason¬ 
able assumption that higher partial-waves are negligible, 
we find 6 '^'^ = (0.45 ± 0.60 ± 0.28)° from the point at 
= 0.14763(51) in [100]i?i. The level at atE^^n = 
0.18712(53) in [ 100 ]il 2 , which is slightly outside the elas¬ 
tic region, yields a phase-shift of = (— 0 . 2 ± 1 . 6 ± 1 . 8 )°. 
It is clear that the P-wave is negligibly small at low en¬ 
ergies. 


We may repeat the analysis described above to deter¬ 
mine the P-wave scattering amplitude, also allowing a 
non-zero P-wave amplitude to influence the spectrum. 
A description of the 22 levels described previously plus 
the [lOOjPi level, using a Breit-Wigner to describe the 
P-wave and a scattering length parameterization for the 
P-wave, k'^ cot S 3 = 1 / 03 , yields 03 = 19(14) x 10® aj, or 
• 03 = 27(20) X 10“®, and P-wave Breit-Wigner pa¬ 
rameters statistically compatible with those given above. 
PT-matrix variations produce similar results, with the P- 
wave amplitude always being compatible with zero. 


B. Coupled-channel tttt, KK scattering 

We now consider the coupled-channel region above 
KK threshold, where tttt — >■ KK is expected to be the 
first significant source of inelasticity. Although we will 
use levels which lie above the 47 r and rjinr thresholds, 
we will not consider those to be open channels. We 
expect the scattering amplitudes featuring these chan¬ 
nels to be very small in the near-theshold energy region 
- experimental support for this assertion comes from 
the measured cross-sections for e+e“ —>■ 47 r [39] and 
e“'"e“ —>■ TyTTTT [40], neither of which has any significant 
value until at least 300 MeV above threshold, likely due 
to the dominance of meson-meson isobars in the ampli¬ 
tudes. Our expectation is that if we were to include 
operators resembling 47 r and/or rjirn into our basis, we 
would extract additional energy levels very close to non¬ 
interacting levels corresponding to weak scattering am¬ 
plitudes, decoupled from the tttt, KK channels that we 
consider. These non-interacting levels will lie somewhat 
above the corresponding thresholds. The formalism to 
understand three-body and higher multiplicity scattering 
amplitudes is not yet complete, although recent progress 
is promising [17, 19, 20]. 

We consider coupled-channel Pl-matrices like those de¬ 
scribed in Ref. [16], using Eq. 4 to define the t-matrix and 
Kij being a 2 x 2 matrix. A particularly useful form for 



where the explicit pole in the first term is an efficient 
way of obtaining a coupled-channel pole in the t-matrix. 
While this parameterization permits a pole to occur in 
the complex energy plane, it is the description of the fi¬ 
nite volume energy levels which determines whether or 
not this pole occurs close to the real axis and is thus 
relevant. We use the Chew-Mandelstam form for the 
phase space, subtracted at the pole position so that 
Re/i(s = m^) = 0. 

We make use of a total of 34 energy levels, shown by 
the black points in Eigures 3 and 4. Eour of these states 
show a significant overlap with a KK operator, whilst 
the remaining levels in the coupled-channel region dom¬ 
inantly overlap with tttt operators. This corresponds to 
using all energy levels below atE^m = 0 . 22 , or below the 
the first unknown “tttt” level, whichever is lowest. This 
spectrum can be described by the Al-matrix of Eq. 6 , 
with = 0 , with parameters 
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m= 0.13170(36)(6) 
g,, = 0.4463(80) (40) 

= 0.71(11)(134) 

= 1.56(94)(30) • ar" 
6.7(26)(143)-ar" 
7 k/?.k/?= 6.8(56)(27)-ar" 


0.20 

-0.24 

-0.27 

0.08 

0 .10' 

1 

-0.77 

0.69 

-0.19 

-0.67 


1 

-0.58 

0.40 

0.90 



1 

-0.03 

0.39 




1 

0.53 


1 


( 7 ) 


We observe a quite reasonable description as measured 
by the X^/A^dof) noting however that some parameters 
are rather strongly correlated, suggesting there is some 
unnecessary freedom in Eq. 6. The phase-shifts, 
inelasticity, g, defined in 


t 


7777,7777 


^KK^KK 

^7nr,KK 


_ I 

gg'^l-^Kk — 1 

“^ipKK 

a/1 - 

2-y/ P-jnrPKK 


( 8 ) 


are presented in Figure 11. We clearly observe the same 
resonant behavior in St^tt in the elastic region that we 
saw previously. We further note that there is very little 
coupling between tttt and KK above KK threshold, and 
that the KK —)■ KK amplitude shows signs of being 
mildly repulsive. That this amplitude describes the finite 
volume spectra rather well can be seen in Figure 12. 

Of course we should be careful not to draw too many 
conclusions from this first description - we cannot be cer¬ 
tain that our choice of parameterization has not forced 
certain features onto the result. To investigate this, 
we consider a range of parameterizations. For exam¬ 
ple we may vary the order of the Tf-matrix polyno¬ 
mial, N, in Eq. 6. We may also consider implement¬ 
ing a “running pole coupling”, where the factors gi in 
Eq. 6 are replaced with energy dependent polynomials, 
Pi (?i(s) = where all g's are real con¬ 

stants. Another variation drops the Chew-Mandelstam 
part of the phase-space, instead just using the simple 
phase-space, Ii{s) = —ipi(s), which satisfies unitarity in 
a minimal way. Such a form is not ideal if we wish to ex¬ 
trapolate far below thresholds, as a kinematic singularity 
appears at s = 0, but we will not have cause to go so far 
below threshold in this case. We summarize these varia¬ 
tions in Table IV showing the resulting y^/Vdof- Fits of 
comparable quality were found representing the elements 
of the inverse AT-matrix as polynomials (as was used in 
Ref [16]), 



^TTTr 


f:KK 


FIG. 11. Coupled P-wave tttt and KK isospin-1 phase shifts 
5 and inelasticity g from a single A-matrix fit. Statistical 
uncertainty shown by the shaded band. The central points 
show the energy levels constraining the amplitude extraction 
with those dominated by TivT-like and gg-like operators shown 
in black and those with significant KK contributions shown 
in orange. On axis circles show the opening of the 47r, KK, 
g-Kn and tvKK thresholds. 


We plot the phase-shifts and inelasticities for a selec¬ 
tion of fits presented in Table IV in Fig. 13 where we see 
that the lattice energy levels very tightly constrain each 
of these forms to give an amplitude description which 
does not vary significantly with parameterization. 


M 

= E (9) 

m —0 

however they produced results with very high degrees of 
parameter correlation, leading to an unreliable estimate 
of statistical error, and we will not discuss them further. 
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Type 


if—matrix 


if—matrix 
withg(s) 


if—matrix 
with 

ii(s) = -ipi{s) 


Explicit form_A^para_X^/A^dof 


11 

_ 9i9j 

- —5- 

—s 

+ 

1 

9kk = 0 

5 

1.75 

Kri- 

9i9j 

—s 

+ 7™ 

5 

'I'KTr ,KK " 0 

5 

1.48 

Kij = 

9i9j 

—s 

1 (0) 
+ Tij 

5 

1kk,kk = 0 

5 

1.37 

II 

9i9j 

—s 

+ 



6 

1.37 

II 

9i9j 

—s 

+I'lf 

s 


6 

1.41 

II 

9i9j 

—s 

1 (0) 
+ Tij 

+ 

g^K = 0 

8 

1.52 

II 

9i9j 

—s 

1 (0) 
+ Tij 

+ 

7™. 

9 

1.39 


K^j 

_ Si(s)9j(s) 1 _,(0). 
rn^—s 'ij ’ 


6 

1.34 


+ 

o' 

II 

1kK,KK = 0, 'yTrn.KK = 0 

K^j 

_ 9 i(!i)gj(s) _ (0). 

m^—s 'ij ’ 


6 

1.33 


+ 

o' 

II 

771-77,7777 = 0, 77777,/C/C ~ ^ 

K^j 

_ 9 i(!i)gj(s) _ (0). 

s ^ij ’ 

9kk{s) = g^^ii 

7 

1.38 


^ \ (0) 1 (1) 
QttttxSj — Qtttt Qmr <5, 

K^j 

_ gi(‘>)gj(s) _(0). 

m^—s fij ’ 

t - „(0) + „(!) s 

' 9kk + 9kk ® 

7 

1.35 



Kij 

_ gi(»)gj(s) ^,(0). 

m'^—s 'ij ’ 


8 

1.37 


/ \ (0) , (1) 
gi{s) =g\ ’ +g\ ’ s 




K.i = 
Kij = 
Kii = 
Kij = 


3iSj 
—s 
9i9j 
ra'^ —s 

9i9j 

rr?,2 —.s 


9i9j 


+ 7^"^- 
^ hj 5 

9kk = 0 

5 

1.57 

+ 7^"^- 
^ hj 5 

^nn ,KK 0 

5 

1.40 

+ 7^"^- 

1kk,kk = 0 

5 

1.58 

1 (0) 

+ Tij 


6 

1.45 


TABLE IV. Coupled-channel if-matrix parameterizations. 




FIG. 12. As Figure 3 with the addition of orange points with 
errorbars showing the spectrum corresponding to the param¬ 
eterization in Eq. 7. 



FIG. 13. As Figure 11 for a range of if-matrix parameteriza¬ 
tions taken from the list given in Table IV. 
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III. RESONANCE INTERPRETATION 

Although we constrain partial-wave t-matrices only for 
real values of energy, either from experimental scatter¬ 
ing, or in this case from finite-volume spectra, the am¬ 
plitudes may be considered to be functions of a complex 
value of s = That the singularity structure of t{s) 
might be important is already apparent if we consider the 
elastic unitarity condition, Imt(s) = p(s) |t(s)p, where 
p{s) = 2fccm(s)/\/s has a square root branch cut begin¬ 
ning at the kinematic threshold. It follows that t{s) also 
has this branch cut and remains single-valued only if we 
consider two Riemann sheets, one where Im k^m is posi¬ 
tive, the “physical” sheet, named because physical scat¬ 
tering corresponds to energies s -I- ie on this sheet, and 
one where Imk^rn is negative, the “unphysical” sheet. As 
more two-body channels open, a greater multiplicity of 
sheets arises, corresponding to the increased number of 
channel momenta. 

The rapid phase and amplitude variation that we as¬ 
sociate with a narrow resonance can be caused by a pole 
at complex values of s = sq = {m — *|r) on unphysical 
sheets®. We may consider our parameterized t-matrices, 
looking for poles at complex values of s, of the form 
tij{s) ^ where we factorize the residue of the pole 
into couplings to each channel, i. 

We find that in every case we considered capable of 
describing the finite-volume spectrum, both in single¬ 
channel and coupled-channels, there is a statistically 
well-determined pole near aty/so = (0.1306— |0.015). 
Parameterizations that do not contain the freedom for a 
resonance pole to occur were not capable of successfully 
describing the finite volume spectra. Figure 14 illustrates 
the position of the found pole, with the lower portion of 
the diagram showing a zoomed region in which the de¬ 
termined pole is shown for a range of different parame¬ 
terizations. A best estimate for the pole position, whose 
uncertainties allow for the spread over parameterizations 
is 




(^0.13055(36) 


1 

2 


0.0150(14) 


decays exponentially as we go lower in energy below the 
KK threshold. 

If we follow the procedure used in previous calcula¬ 
tions, making use of the n baryon mass determined on 
these lattice configurations, to set a physical scale we find 
at = ((( , where at mo is determined using lattice QCD 

computation and = 1672.5 MeV is the experimen¬ 

tal mass. Using 16 distillation vectors on this lattice we 
have determined Otmo = 0.2789(16), which leads to an 
approximate pion mass of = 236 MeV. 

With this scale setting, the resonance pole is located 
at y/'^ = [783(2) — |90(8)] MeV. The scale-set Breit- 
Wigner mass and width of Eq. 3 are mew = 790(2) MeV, 
Fbw = 87(2) MeV, and a plot of the corresponding 
phase-shift with the scale-set energy is presented in Fig¬ 
ure 15. 

An earlier calculation by the Hadron Spectrum Col¬ 
laboration considered elastic tttt scattering using lattice 
configurations with the same quark and gluon action, but 
with a larger value of the u, d quark mass, such that 
the pion had a mass of 391 MeV [22]. We compare the 
Breit-Wigner parameters in this study with those deter¬ 
mined for m^n- = 391 MeV in Table III^. The correspond- 


0.08 



- 0.01 



-0.04 




- 0.012 


-0.013 


-0.014 


-0.015 


-0.016 


-0.017 




0.128 0.129 0.130 0.131 0.132 0.133 



The corresponding coupling to the tttt channel also shows 
very little variation under parameterizations with a good 
estimate being a* = 0.049(3) The cou¬ 

pling to KK, which only arises in coupled-channel anal¬ 
ysis is not well determined, having a large statistical un¬ 
certainty. Along with the observation that the elastic 
data can be very well described without invoking any 
KK amplitude, we conclude that we have not reliably 
constrained . This is to be expected as the effect 
of the KK amplitude on the spectrum in finite-volume 


-0.018 


FIG. 14. Resonance pole position on first unphysical sheet. 
Zoomed region shows the pole found for a variety of pa¬ 
rameterizations: various descriptions of the elastic amplitude 
(red), various A-matrix descriptions of the coupled-channel 
t-matrix, using the Chew-Mandelstam phase-space (orange), 
and using the simple phase-space (blue). 


a conjugate pole must also be present at Sq, but this pole is 
usually much further from the region of physical scattering. 


^ The results presented in [22] suffer from a small error in the 
computation of the off-diagonal data covariance, which we fix 
here, leading to a very small shift (at the level of Itr) in the 
quoted Breit-Wigner parameters with respect to that reference. 
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1 vn 



FIG. 15. Elastic Breit-Wigner fit of Figure 8 with the lattice 
spacing set as described in the text. 


This work Ref. [22] 



0.03928(18) 

0.06906(13) 

atrriR 

0.13175(35)(5) 

0.15095(14)(4) 

9R 

5.688(70) (26) 

5.698(97)(3) 

m,r 

236(2) MeV 

391(1) MeV 

rriR 

790(2) MeV 

855(1) MeV 

Re(^^) 

783(2) MeV 

853(2) MeV 

-2Im(y^) 

85(2) MeV 

12.4(6) MeV 

1 Cttti- I 

288(4) MeV 

162(4) MeV 

Arg(c,r7r) 

-0.059(1) TT 

-0.033(1) TT 

X^Ndot 

^ = 1-25 

IS = 0.98 



FIG. 16. Elastic P-wave tttt isospin=l phase shifts for pions 
of mass ~ 391 MeV (blue) and ~ 236 MeV (red). Gurves 
correspond to two parameter Breit-Wigner fits as described 
in the text and summarized in Table III. 


ing phase-shifts are shown in Fig. 16. We note that the 
Breit-Wigner couplings gji show good agreement between 
the two different quark masses, as has been suggested 
in unitarized versions of chiral perturbation theory [42]. 
The apparent difference in the pole residue coupling, 
is completely explicable in terms of the P-wave barrier 
- since the t-matrix near threshold must behave like 
t{s) ^ we may consider k{so) where k{so) is 

the cm-frame momentum at the pole position. It follows 
that 


c 


391 

TTTT 


^236 

'-'TTTT 


0.56(2) 


^236 

173 - 7.7i 
312 - 26.6f 

0.552, 


and thus 


=2361 


TABLE V. A comparison of the results of this study and 
ref. [22]. These numbers compare the Breit-Wigner descrip¬ 
tion only and the quoted pole is from that single parameteri¬ 
zation. 


IV. SUMMARY 


The p resonance was extracted using a detailed spec¬ 
trum of lattice QCD energy levels working in a (~ 4fm)^ 
cubic volume with a pion mass of 236 MeV. Using the 
variational method and a large diverse basis of operators, 
energy levels were obtained in the elastic tttt scattering 
region and the near-threshold coupled-channel tttt — KK 
energy region and these were used to constrain the 7=1 
= 1“ and = 3“ scattering amplitudes. By 
making use of the formalism relating the elastic and 
coupled-channel scattering amplitudes to the spectrum 
of eigenstates in a finite volume, we were able to extract 
phase-shifts and inelasticity for the coupled tttt — KK 
system. The elastic region was found to feature a nar¬ 
row resonance, which persists when the coupled KK 
channel is also considered. A range of t-matrix param- 
eterizations lead to consistent resonance parameters in 
the sense of a pole in the complex energy plane. Us¬ 
ing the il baryon to set the scale this pole is located at 
= (783(2) — |90(8)) MeV. A simple Breit-Wigner 
description of the elastic amplitude works well and gives 
a coupling that is consistent with the value determined at 
a larger pion mass and that extracted from experimental 
data. 

In a coupled-channel analysis we found the tttt — KK 
system to be only weakly coupled for the range of energies 
we considered, and only small phase-shifts were observed 
in the KK channel. This is only the second example 
of the extraction of a coupled-channel scattering matrix 
from lattice QCD, following the earlier study of ttK, rjK 
[15, 16]. This first exploratory study above KK threshold 
neglected three- and four-hadron contributions that have 
been observed to be suppressed in experimental studies. 
Progress in the development of a finite-volume formalism 
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capable of dealing with these higher-multiplicity channels 
is ongoing [17, 19, 20]. 

The consistency and broad applicability of these meth¬ 
ods to extract resonance properties from lattice QCD is 
now being demonstrated with successful applications at 
multiple pion masses in various quantum numbers. Fu¬ 
ture studies will aim to shed light on longstanding mys¬ 
teries such as the ao(980), /o(980) resonances, and inves¬ 
tigate excited hadron states with exotic quantum num¬ 
bers. 
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Appendix A: Scattering analyses with an incomplete 
spectrum 


In our basis we included operators that were specifi¬ 
cally constructed to resemble both the tttt and KK states 
one would expect to exist in the absence of meson-meson 
interactions, as well as fermion bilinears, which re¬ 

semble qq-like constructions. As is visible in Figure 6, 
the low-lying states overlap with both sets of operators. 
This is what we might expect for a system containing 
a resonance that is dominantly qq, but coupled to the 
decay channel tttt. 


In order to gauge the importance of including TTTr-like 
operators in the basis, we perform variational determina¬ 
tion of the spectrum using only the “single-meson-like” 
operators in each irrep. The energies so determined are 
converted to elastic phase-shifts and are plotted in Fig¬ 
ure 17. For comparison, we also show the Breit-Wigner 
elastic phase-shift curve determined in the text. As 
we have previously suggested [45], when a narrow res¬ 
onance is present, using only “single-hadron-like” oper¬ 
ators tends to provide energies which lie roughly within 
one hadronic width of the mass of the state. It is clear 
that this limited operator basis may be useful to suggest 
the presence of a narrow state, but it cannot determine 
the resonant properties of such a state. 



FIG. 17. Elastic phase-shift points (colors and symbols as in 
Figure 7) corresponding to the spectrum extracted using only 
“single-hadron-like” operators, tpT^. The grey band shows 
the elastic Breit-Wigner amplitude determined from the full 
spectrum as described in the body of the manuscript. The 
blue dashed vertical line and band indicate the Breit-Wigner 
mass and width. 






14 


Appendix B: Flavor structure of KK operators 


Our “single-meson-like” and isospin-1 “TTTr-like” oper¬ 
ators have already been described in detail elsewhere 
[21, 28-30]. The isospin-1 “ATAT-like” operators are con¬ 
structed to have positive G-parity (corresponding to neu¬ 
tral states with negative G-parity). For example, the 


Iz = -1-1 component is proportional to, 

(sFiu) {dT2s) + G [(sFiw) {dT2s)] , 

where encodes the spin, derivative and momentum 
structure of the operator. Here G is the G-parity trans¬ 
formation, G [ 9 irg 2 ] = G G[( 72 ]rG[gi], where G is the 
G-parity of the underlying spin and derivative structure, 
i.e. that for a flavorless qTq operator. The KK operators 
are projected onto definite irreps of the relevant symme¬ 
try group by summing over relative momenta in exactly 
the same way as for tttt operators. 
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